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Abstract. The discovery that the Galactic centre emits flares at various wavelengths represents a puzzle concerning their 
origin, but at the same time it is a relevant opportunity to investigate the environment of the nearest super-massive black 
hole. In this paper we shall review some of our recent results concerning the tidal evolution of the orbits of low mass satellites 
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phenomena Uke near infra-red and X-ray flares from Sgr A*. 
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INTRODUCTION 

It is nowadays widely accepted that most galaxies, if not 
all, contain a super-massive black hole at their centre. 
How these black holes are formed and how they are re- 
lated to galaxy formation is still an open issue, see e.g. 
[1]. The nearest super-massive black hole (SMBH) is 
located at the centre of our Galaxy, Sgr A*. Stellar or- 
bits determinations [2] show that a central dark mass of 
(3.7±0.2) X lo6[/?o/(8kpc)]3M0 (with an updated value 
of (3.61 ±0.32) X IO^Mq from observations with SIN- 
FONI [4]) is confined within a radius of 45 AU. The 
proximity of the Galactic centre (« 8 kpc) provides us 
a unique opportunity to study in detail the environment 
of massive black holes (see e.g. [5] that contains an ex- 
haustive description of the Galactic centre), and possi- 
bly of their formation. Due to its proximity, Sgr A* is 
the only galactic nucleus that can be investigated using 
VLBI on sub-AU Unear scales, and observations in 2007 
have confirmed structure in Sgr A* on scales of just a 
few Schwarzschild radii [6]. 

The recent detection of rapid flare activity coming 
from Sgr A* represents additional strong evidence for 
the presence of a massive black hole. Their time scales 
strongly suggest motion only a few Schwarzschild radii 
away from the central black hole. If such flares are un- 
correlated single events of accretion, we believe that the 
most plausible candidate for their energy source is the 
tidal part of the gravitational potential energy. We sug- 
gest that Galactic flares are produced by the final ac- 
cretion of a single dense object, like an asteroid or a 
comet, with a mass of lO^^g. We explore this sce- 
nario, and we show that the light curve of a flare can 
be deduced from dynamical properties of geodesic orbits 
around black holes and that it only weakly depends on 
the physical properties of the source. 



FLARES FROM THE GALACTIC 
CENTRE 

Flares from Sgr A* have been detected in X-rays (Chan- 
dra and XMM-Newton satelUtes) and infra-red (VLT 
adaptive optics imager NACO, SINFONI infrared adap- 
tive optics integral spectrometer on the ESO VLT). The 
description of this flaring activity can be found in several 
papers, e.g. [7, 8, 9, 10, 11, 12]. These flares are modu- 
lated on a short time-scale with average periods of w 20 
minutes, and multi-wavelength campaigns found that the 
time lag between X-ray and NIR flare emission is very 
small, strongly suggesting a cormnon physical origin, see 
e.g. [13, 14, 15, 16]. 

The main features of the observed flares are: a) du- 
ration: fa l/4h - 2h; b) frequency: few events per day; 
c) energy release: of the order of 10^^ erg/sec; d) linear 
polarization (NIR, radio): of the order of few to few ten 
percent, with stable polarization angle; e) strong quasi- 
periodic modulation on time scales of minutes. 

Several different models were proposed to explain 
these phenomena: disk instabiUties (e.g. [17, 18, 19, 20]), 
star-disk interaction (e.g. [21]), expanding hot blobs (e.g. 
[15, 22]), hot spot/ring models (e.g. [16, 13, 14, 10, 
12, 23, 24, 25]). Here we do not want to compare and 
contrast the above proposed model. Our aim is to propose 
a scenario in which Galactic flares are produced by the 
final accretion of a single dense object, like an asteroid 
or a comet. The main motivations of this approach are: 
a) the time-scales of flares strongly suggest motion only 
a few Schwarzschild radii away from the black hole; b) 
the energy release during these flares corresponds to the 
source mass of the order of lO'^'^g, see e.g. [8]; c) it seems 
reasonable to assume that stars at the Galactic centre 
are surrounded by planets and by other small orbiting 
bodies, like asteroids and comets (hereafter: LMS, low- 



mass satellites); d) tidal evolution of the orbits of LMS 
will lead to capturing orbits; e) tidal effects in the vicinity 
of a black hole can release a significant amount of energy. 

Asteroid-like objects fit both observational facts (time- 
scales and energetics): they are of the right mass and they 
are tidally disrupted closer to the black hole than gaseous 
blobs. In active galaxies, such events would be outshined, 
but in the inactive centre of our Galaxy, it is possible to 
observe them. 

In the following Sections we shall expand on the above 
motivations. Most of this work is based on our recent 
papers, e.g. [26, 27, 28, 29], where a thorough discussion 
and more details can be found. 



LMS AT THE GALACTIC CENTRE 

According to [30], the Edgeworth-Kuiper belt of our So- 
lar System may still (4.5 billion years after its creation) 
contain as many as 2 x 10*^ objects of radii < 10 km. 
It appears reasonable to assume that stars at the Galac- 
tic centre are surrounded by planets and by other small 
orbiting bodies (LMS), like asteroids and comets, that 
may be stripped off their parent stars by tidal interaction, 
while approaching the black hole. Therefore, we expect 
that there must be a considerable number of solid objects 
that cluster the Galactic centre. 

Little is known about dynamics that determines the 
fate of such low mass satellites. However, some hints 
come from the investigation of core dynamics at the 
Galactic centre [32]; from the investigation of mecha- 
nisms that may cause a star to move inwards into a mas- 
sive black hole, in connection with possible detection of 
gravitational waves with LISA, see [33] and references 
therein; from detailed dynamical studies of the evolution 
of the Solar System, see e.g. [34, 35, 36]. 

One might expect to find stellar system satellites all 
the way down to the black hole at the Galactic centre, 
with a fair proportion of them on low-periastron, highly 
eccentric orbits. A significant work will be done by tidal 
forces near the low periastron. This work partially dis- 
sipates as heat, and as a result it lowers orbital energy 
and may start significant tidal evolution of the orbit, see 
[27]. When periastron crossing frequency and fundamen- 
tal quadrupole frequency of the object are the same, that 
is at resonance, tidal interaction is strongest. In this con- 
text, one must consider two classes of LMS: those that 
are gravity dominated and those that are sohd state forces 
dominated. 

For gravity dominated objects the fundamental fre- 
quency is Vg ~ l^jGp /3n, while for solid state domi- 
nated objects it is Vs = \cs/R (p is the density of the 
body, R its radius and Cj the speed of sound). Taking 
as typical values Cg ~ 5 km/s and p ~ 5 g cm"^, we 



find that the radius dividing the two classes is about 
1000 km, i.e. the radius of the asteroid Ceres. Accord- 
ingly, all gravity dominated satelhtes must have about 
the same fundamental quadrupole frequency, that cor- 
responds to a period of about 1 hour and all smaller 
satellites have shorter fundamental periods. Thus, gravity 
dominated satellites start a rapid, resonant tidal evolution 
when their periastron reaches rp w 10 rg (rg = GM^h/c^), 
while solid state dominated bodies may undergo signifi- 
cant tidal evolution even closer to the black hole. 

The difference between these two classes of objects is 
better understood introducing the Roche radius and the 
effective Roche radius. The effective Roche radius is the 
radius at which cOqT} = 1, where cOq is the fundamen- 
tal resonant frequency of quadrupole modes and 7f = 
(2rp/GMbh)'''^ is the periastron crossing time [26]. In 
Fig. I we show the Roche radius and the effective Roche 
radius as a function of object size, mass and density for 
different bodies expected to populate the Galactic cen- 
tre. This graph accentuates the very particular position 
of asteroids. With the exception of compact stars such as 
white dwarfs, neutron stars or stellar mass black holes, 
they are the only objects that can survive the tidal field of 
the black hole when held together by soUd state forces. 
All other objects such as planets, main sequence stars and 
of course giants find themselves beyond the Roche radius 
above ten or a few ten gravitational radii. Note, however, 
that molten asteroids become gravity dominated. If melt- 
ing occurs suddenly below ~ 10 rg, the phase change 
results in a dramatic rapidly developing tidal disruption. 
We intend to show that why and how this may happen. 

TIDAL EVOLUTION OF LIVIS ORBITS 

Tidal effects on a star encountering a SMBH have been 
investigated by several authors, see e.g. [26, 31] and 
references therein. However, tidal forces will not affect 
only the shape and the structure of the in-falling body, 
but are responsible for the evolution of the orbits and the 
rotation periods of the involved objects. 

The first stages of tidal evolution of LMS orbits can 
be investigated using Hut's formalism described in [37, 
38, 39]. This formalism is suitable for describing slow 
tides, i.e. tides that change slowly with respect to the 
frequency of the fundamental quadrupole modes of the 
object that tides are acting upon. In this case the relevant 
forces that exchange angular momentum and energy can 
accurately be calculated by approximating the deformed 
object adding two additional bulge masses pL on the 
surface of the star. In this case tidal evolution of orbits 
depends only on the parameter a = Iq/sq, which is the 
ratio of orbital angular momentum of the system and spin 
of the secondary at equilibrium. Therefore, completely 
different binary systems with the same value of a behave 
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FIGURE 1. Astrophysical objects capable of releasing sig- 
nificant tidal energy before falling down the Galactic centre 
black hole. The positions of different celestial objects are also 
marked in this diagram. Note that the effective Roche radius for 
solid state dominated asteroids is below 2 rg, while the Roche 
radius for these objects is about 10 rg. 



in the same way. See [28] for details. 

We find that the evolution of orbits strongly depends 
on the initial value of ip o, which is the ratio between 
the periastron distance ip and the radius of the equilib- 
rium circular orbit ao. If an orbit starts with Fp o below 
some critical value, it will always end up in in-spiralling, 
regardless of the value of a. Orbits starting above this 
critical value of rp o will either circularize, or become un- 
bound. 

Tides induce forced oscillations on the object, excit- 
ing its eigenmodes. They are non-resonant, as long as 
the frequencies co of spectral components of tidal inter- 
action are much lower than the quadrupole eigenfrequen- 
cies (Oq of the object. If, on the other hand, spectral com- 
ponents of tides are strong at frequencies close to eigen- 
frequencies of the quadrupole oscillations of the object, 
which happens when the orbit touches the Roche radius, 
tides become resonant. Fig. 2 and Fig. 3 summarize our 
results. 

As indicated in Fig. 2, tidal evolution of the orbit has 
two distinct phases: the non-resonant phase, well de- 
scribable with Hut's slow tides formalism and the res- 
onant phase, which was studied numerically in affine ap- 
proximation by [3] and in Newtonian approximation by 
[26]. The energy loss per periastron passage can in both 



FIGURE 2. As tides drive the tidal bulge across the surface 
of the solid asteroid, work deforming the body is done at the 
expense of orbital energy; little angular momentum is deposited 
in spinning the body because of its low moment of inertia (non- 
resonant phase). When the body melts and spends more and 
more time beyond the Roche radius, it elongates, increases the 
moment of inertia which allows both orbital angular momen- 
tum as well as orbital energy to be transferred to spin (resonant 
phase). 



cases be described by the expression 

where e^{j5eff) is the eccentricity of the ellipsoid to 
which the originally spherical object is deformed by 
tides; it depends essentially only on the effective Roche 
penetration parameter jSg// = r^^yy/rp'. In the case of 
gravity dominated objects jSg// = j3 = rR/rp and in the 
solid case jSf.// is smaller; it can be expressed as jSg// = 

m\[%7p ~ 3300 ^ ^h^i'^ 3300 

comes from inserting the approximate value for the ve- 
locity of sound and the density of solid bodies. The 
time scales for non-resonant and resonant tidal evolution 
phase can then be calculated (cf. [27]). They can be ex- 



This approximation is valid for small amplitude resonant tides and is 
valid as long as most part of the orbit is outside the Roche radius, i.e. on 
highly eccentric orbits whose periastron reaches below Roche radius. 
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FIGURE 3. Resonant and non-resonant time-scales as a 
function of j8 for p = po, R = lOKm, Q = 200. j3 is tiie Roche 
penetration parameter fi = rn/rp and Q is tfie resonant damping 
factor. 



pressed in the form 

2 - 

10 fER{e) / 10km\ /po^« 



fER = 2.9 X lO'V X 



(3) 



where Q is the resonant damping factor, fEN{e) and 
fER(e) are factors of order unity depending on the ec- 
centricity of the orbit (e), and e^(j3) is a function that 
peaks at resonance i.e. at j3 = 1. 

These time-scales are plotted in Fig. 3 for 10 km ob- 
jects with p ~5g cm^^ and 2=200. Also plotted in this 
figure is the approximate melting time, i.e. the time it 
takes to accumulate enough tidal work inside the object 
to melt it {tneit = E,„ei,t,„hl^\ E,„ei, ^ 10" '°m c^), and 
the thermal diffusion time for such an object (assuming 
the thermal diffusion coefficient D ^ 10"''m^5"^ char- 
acteristic of most rock). It is clear from this figure that 
10 km objects melt for j3 > 1, i.e. if their periastra reach 
below 33 rg in less then 3 x 10^ years. Asteroids that 
do not penetrate as deep or are smaller would not melt, 
because their thermal diffusion time scale is too short^. 
On the other hand asteroids that are scattered or released 



Note that tidal heating is rather slow, which indicates low tidal stress, 
so that the break-up of such asteroids is not expected. 



to highly eccentric orbits with pariastra at 10 rg melt in 
only a few years. When melted, asteroids become reso- 
nant and their tidal evolution time-scale shortens consid- 
erably, certainly to below a few billion years but possibly 
even to a much shorter time scales, depending on their 
size and most of all on the periastron of the orbit that 
they were ejected or scattered to. Increased tidal heating 
does not essentially raise their temperature, because con- 
vection sets in so that low temperature (~1000 K) black 
body radiation can dissipate tidal heat. One can conclude 
that all asteroids found within 33 rg are melted, non- 
evaporated and are bound to end up in the black hole. 



THE FREE FALL PHASE 

At present we do not have a complete theory describing 
all the steps of evolution toward plunging into the black 
hole. Clearly, at 33 rg asteroids have entered the resonant 
part of their tidal evolution in relativistic regime. Fig. 2, 
based on Hut's slow tide approximation, indicates that 
the characteristic of this part of evolution is the trans- 
fer of orbital angular momentum to spin. In the above 
resonance regime this transfer can not be effective by 
spinning the object faster, but by deforming it into a 
longer and longer thread, thus increasing its moment of 
inertia. Some understanding of this mechanism can be 
obtained considering the viiial theorem applied on the 
self-gravitating mass of an orbiting asteroid. In [26], the 
authors consider a self gravitating object moving in the 
gravitational field of a point mass (Newtonian approxi- 
mation). Following the steps similar to those in deriving 
the virial theorem, they obtain the equation 
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where Wi„t is the internal energy of the asteroid, is 
the self-gravitational potential energy, R is the orbital 
vector from the black hole to the asteroid, Q is the 
quadrupole tensor of the deformed asteroid, and J is the 
momentum scalar (J = J pr^dV). If one approximates 
the deformed asteroid by an axis-symmetric ellipsoid 
with the symmetry axis in the direction of «, then Q can 
be written in the form Q — q{3nn — I), where I is the 
identity matrix and q = — |M* ( rj^^^j^ ~>'1g)- Furthermore, 
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Neglecting the equatorial radius with respect to the polar 
one {rpoie r^q), this becomes 

1 M* ( J J , , ^ 1 9 \ 

"^int + = ^ ( C0|r^„;,(6cos" a - 2) + J , 

(6) 

where a is the angle between R and h and 0)^: = 
Gmb\i/R^ is the Keplerian angular velocity correspond- 
ing to a circular orbit at radius R. In stationary equilib- 
rium the right hand side of this equation vanishes and the 
expression reduces to the well known virial theorem. In 
the other extreme, if tides had elongated the object a few 
times its original radius, the left hand side of the above 
equation becomes negligible with respect to the first term 
on the right, so that the equation has become a second 
order linear differential equation for Its solutions 

are oscillatory if the average value of (6cos^ a — 2) is 
positive, and exponentially expanding if it is negative. 
For slow tides, the bulge points almost in the direction 
of R (a ^ 0) and accordingly r^^^^ is an oscillatory 
fiinction. However, above resonance, the tidal bulge may 
point more and more in the direction of motion, which 
changes the character of the above equation making the 
tidal bulge grow exponentially. 

The virial theorem is not sufficient to describe this 
phenomenon in any detail, because it does not predict the 
value of the angle a. Therefore, we study the last stages 
numerically with a fully relativistic code. This code, de- 
scribed in [29] builds an image of the tidally evolving 
object assuming that it starts as a molten sphere. In calcu- 
lating the developing shape of the body fully relativistic 
dynamics is taken into account, with only one simplifi- 
cation that the intemal pressure can be neglected with 
respect to the high kinetic energy density of flow in the 
tidal bulge. Evolving images of the object are constructed 
as they would be seen momentarily in the frame of a dis- 
tant observer. A selection of a few slides in Fig. 4 shows 
such a development and confirms predictions of expo- 
nential growth in the above simple argument. In Fig. 5 
we also show the development of the visible length of 
an originally spherical object as it is stretched by tides, 
obtained in a number of numerical experiments. We find 
that resonant tides develop exponentially growing tidal 
tails which are fastest and most pronoimced in the vicin- 
ity of the black hole. 

At this point it is important to appreciate the difference 
between different tidal evolution problems. In the Kep- 
lerian case, described so well by Hut's mechanism, the 
tidal dissipation of orbital energy and exchange of an- 
gular momentum between orbits and spins brings orbital 
parameters closer and closer to equilibrium at the min- 
imum of the effective potential. This equilibriiun may 
be reached (depending on initial conditions) either at a 
distance where the pair is well separated so that no part- 
ner is within each others Roche radius, or on an orbit 



where one partner finds itself beyond the Roche radius 
of the mutual potential. In the first case, tidal evolution 
has stopped since tides can no more dissipate energy (the 
only dissipation mechanism left is by gravitational radi- 
ation which is much slower in most cases). In the other 
case, the member that finds itself beyond its Roche ra- 
dius spills over to the companion by building an accre- 
tion disk around it. Unless there is a strong interaction 
with other bodies in the vicinity, the accretion phase de- 
velops, according to Hut's mechanism, during a stage 
when the members of the binary are (almost) co-rotating. 
Therefore, according to the above discussion, the angle 
a remains close to zero, keeping equation 6 in oscilla- 
tory regime. This ensures that the shape of the mass- 
losing star remains stable. Mass flow through the iimer 
Lagrangian point is made possible only by the (slow) 
transfer of angular momentum enabled by energy dissi- 
pation in the accretion disk. If, on the other hand, the 
accreting mass is a black hole and the inner Lagrangian 
point is close to the maximum of the black hole effective 
potential, then most geodesic orbits from the mass-losing 
companion lead directly behind the horizon of the black 
hole. No stationary accretion disk can form in this case 
and the mass-losing companion is sucked to the black 
hole in one gulp as a characteristic exponentially elongat- 
ing string. In our numerical experiments, we found that 
the dividing line between non-relativistic and relativistic 
accretion is well represented by a relativistic parameter 

c 

_ E/mc^ - Vmin 
^ max * mm 

where E is the orbital energy of the asteroid and Vmax 
and Vmin are the maximum and minimum of the rela- 
tivistic effective potential. Relativistic accretion occurs 
only if ~ 1 . Even orbits winding beyond the Roche ra- 
dius and close to the black hole, but with orbital energy 
near to Vmin, do not accrete relativistically. The elonga- 
tion of the object does not develop exponentially but lin- 
early, as shown by the curve 'a' in Fig. 5. The meaning 
of ^ can be visualized in Fig. 6. It is clear from this fig- 
ure that objects on bound orbits accrete relativistically 
only if their orbital angular momentum decreases below 
AmrgC. Thus, relativistic accretion can occur only from 
orbits with periastra between 4 and 6 Vg. Such unstable 
relativistic orbits wind about the black hole few times 
and unwind from there into the black hole. The winding 

periods {T = /rg) on such orbits are between 

15 and 28 minutes (for a 3.6 x lO^M© GC black hole). 

The main conclusion of this section is the realization 
that tidal evolution of binary systems with a black hole 
has three, and not only two possible outcomes as in 
Newtonian theory. Namely, accretion of a smaller object 
onto a black hole would not occur by gradual Roche 
lobe overflow as in a contact binary, but in a single 



FIGURE 4. A selection of images from simulations of tidal development of an 0.03 rg object following an elliptic orbit around 
a black hole observed from above the orbital plane (frames on the left) and as seen from 1° above the orbital plane. Captions to 
left of frames indicate the size of grid spacings; green grid lines mark coordinates which are integers of rg. Effects of gravitational 
lensing are clearly seen in frames on the right. 




relativistic gulp, if the Roche lobe overflow has to occur 
at the inner Lagrangian point between 4 and 6rg from 
the black hole. The energy released during such an event 
is probably comparable or more then the energy released 
during the whole disk accretion phase and the light-curve 
is modulated with the last period of revolution about the 
black hole. The new relativistic outcome occurs, if the 
binary started as detached and the Roche lobe overflow 
occurs suddenly (due to melting) near the periastron 
which is below 6r„. 
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FIGURE 5. The length 

^^pole of the image of a tidally- 
evolving object seen above the orbital plane as a function of 
time for different orbital parameters: a) elliptic orbit close to 
innermost stable circular orbit (ISCO) with = 6.2122 rg, 
rp = 6.2097 rg and ^ = 2.44 x 10^^; b) plunging parabolic 
orbit with E/mc^ = 1, l/mrgC = 3.999998 and = 1 + 10"*; 
c) plunging hyperbolic orbit with r(V,^ax) = 3.6 rg and ^ = 
1 + 10"*; d) elliptic orbit close to ISCO with = 6. 1 rg, rp = 
5.95 Tg and f = 1 - 3 X 10"*; e) elliptic orbit with = 20 rg , 
rp = 4.45 rg and = 1 — 3 x 10"*; f) scattering parabolic orbit 
with £/mc2 = 1, l/mrgC = 4.000002 and = 1 - 3 x 10"*; g) 
plunging hyperbolic orbit with E/mc^ = 3.0 and l/mrgC = 6 
and C, = 6.35; h) linearly growing deformation. In cases b,c,e 
and f, the simulation was halted before the object evolved 
beyond a limit of numerical accuracy. 
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FIGURE 6. The Schwarzschild (black) and Kepler (light 
gray) effective potentials for different values of dimensionless 
angular momentum (A = ^^ff-^) 



THE ROLE OF MAGNETIC FIELDS 

NIR flaring emission from Sgr A* is believed to be pro- 
duced by synchrotron radiation of electrons with Lorentz 
factors 7g ~ 10^ — 10^ moving in a magnetic field of 
10- 100 G [41, 7, 42]. In this section we would like to 
show that the development of such strong magnetic fields 
may be a natural development during the short event of 
relativistic accretion. We would also like to show, in a 
qualitative manner, that if certain physical conditions are 
met, then the global magnetic (and electric) field may 
play the decisive role in changing the large kinetic energy 
of the exponentially extending object into kinetic energy 
of electrons and finally into radiation, released just be- 
fore the object forever disappears behind the horizon of 
the black hole. 

There is still no general agreement on the strength 
of the magnetic field at the Galactic centre. In [43] it 
is pointed out that the field may be strong (of order 
1 mG) and globally organized (see e.g. [44]) or globally 
weak (of order 10 flG) with regions where it is stronger 
(see e.g. [45, 46]). In any case, the observed strength is 
orders of magnitude lower than requested by synchrotron 
models. Assuming that the object is permeated with such 
low magnetic field, our numerical simulations suggest a 
simple mechanism to increase it, as well as to provide 
electrons with high y^- 

We have found that objects moving along ^ ^ 1 or- 
bits end up being exponentially squeezed and elongated. 
Since the inward squeezing velocity is low (v/c ^ 10"^), 
it is possible to use the magneto-hydrodynamic (MHD) 
approximation to describe the evolution of the magnetic 
field 
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Vx (VxB) = -— + Vx (v?xB) 
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(8) 



where a is the electrical conductivity and v is the 
plasma velocity. Assuming typical values a = 4.14 x 
10* n"'m"' for molten aluminium and L = 100 km for 
the size of the object, the magnetic Reynold's number 
becomes very large 



R,„ = lioavLKi 10^ 



Therefore we may neglect the diffusion term in the above 
equation 

dB 

^ = Vx(vxB), (9) 

which, according to Alfven's frozen-flux theorem, leads 
to magnetic flux conservation. Consider an object with 
initial cross section which is permeated with the ini- 
tial magnetic field Bq- If such an object is exponentially 
squeezed, magnetic flux conservation increases the mag- 
netic field density as B = BqSo /S. Since the tidal defor- 
mation tensor is traceless, the average cross section area 
of a squeezed and elongated object is inversely propor- 
tional to its length (in the GR case this result is only ap- 
proximate). Therefore, the magnetic field at a later times 
can be estimated as 



rpolejt) 



where Rq is the initial radius of the object and r^o/g (t) the 
major axis of the ellipsoid, to which the object has been 
deformed at time t. In Fig. 5, results of dynamical sim- 
ulations show length increase to follow the exponential 
law up to five orders of magnitude. The limit comes from 
numerical limitations of our simulation only - at this fac- 
tor of expansion the numerical grid became too coarse 
to trust further calculations. In fact, in the approximation 
taken, where internal pressure was neglected, (almost) no 
limit to the expansion factor is expected. Equation 10 ob- 
viously indicates that neglecting he internal pressure has 
its limits. When the magnetic pressure becomes compa- 
rable to dynamical pressure involved in driving the ex- 
panding flow, the increase of magnetic pressure must be 
stopped, impeding further contraction. 

In order to estimate dynamical pressure inside the 
exponentially expanding object, we consider a simpli- 
fied problem with a solution that is similar to the expo- 
nentially expanding asteroid solution from previous sec- 
tion. Consider a spherical self-gravitating incompress- 
ible fluid that suddenly finds itself in a strong field of 
a cilindrically symmetric gravitational quadrupole ((^^ = 

2 

'^{x^+y^ — 2z^)). An exact solution of this hydrody- 
namic problem can be found such that the sphere is con- 
tinually deformed into and exponentially extending el- 
lipsoid with axes rpoie{t) = Rqc'I'^ and req(t) = Rqc^'^'^'^, 
with T = (2rg/r^)~^/-^/c, and the central pressure de- 
creases (for t 9 rpoie '> req) as 
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(11) 



Note that the central pressure is 3/4 of what it would be 
in hydrostatic equilibrium. If magnetic pressure pushing 



outward at equator is added and balances the central pres- 
sure, there is no more pressure gradient toward equator, 
so radial compression must stop. The negative magnetic 
pressure in the axial direction is also working against ex- 
pansion, however the expansion velocity along this axis 
is much larger than the compression velocity, so expan- 
sion can only be halted if the kinetic energy of expan- 
sion is transformed into some other form, as for example 
magnetic energy or kinetic energy of electrons. Equat- 
ing the central hydrodynamic pressure and the magnetic 
pressure = 2^^^' '^^^ obtains the equation for the 
equatorial radius when compression stops {Req) 
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Inserting the numbers used in previous arguments (T 



11 200sec, p = 5g cm ^) and measure Rq in units of 
10 km and Bq in units of mG, then the above becomes 
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(13) 



and the ratio of polar (Rpoie) to initial radius follows as 



R 
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Ro V^o[10km]y 



1/3 



(14) 



Equations 13 and 14 tell that a 10 km asteroid permeated 
with an initial magnetic field of 10 mG stops compress- 
ing when its equatorial radius is 27 meters, the magnetic 
field has increased to 133G, the polar radius extended 
to 1.33 X lO^km, and the whole event lasted ^ AOmin - 
less than 3 turns around the black hole. The final veloc- 
ity of the pole with respect to the centre of mass is about 
6.6 X lO^km/sec and the kinetic energy of expansion of 
the order 2 x lO^^m c^. 

These 40 minutes are rather eventful. Very few pro- 
cesses in the universe release so much energy so fast. 
The work of gravity on the still extending object can go 
on to extract in a few more turns up to ^ O.lmc^ in tidal 
energy ([26]). 

This mechanism works as long as the magnetic dis- 
sipation time-scale Tq = ajJ^r^^ is longer than the dy- 
namical time-scale. Initially Td/sjw1.2 x lO^sec and this 
condition is certainly well satisfied. Finally, when the ob- 
ject is only 27 meters in diameter, the value of Tq drops 
to ~ 15min if the same molten metal conductivity is as- 
sumed. It is quite clear that this assumption can no longer 
be respected, since after the release of so much energy, 
the object can no longer be considered a molten metal 
but has obviously turned into hot plasma. The conductiv- 
ity probably increases and keeps magnetic flux constant. 

As long as the magnetic field is low, it foUows the 
motion of the plasma. However, when it reaches high 



enough values, matter starts to move according to this 
magnetic field, thus providing a mechanism for acceler- 
ating electrons to relativistic velocities. The gravitational 
force which guides the motion of positive charges is neg- 
ligible for electrons. They are pushed out by magnetic 
pressure until charge separation induces enough electric 
field to keep them bound to positive charges beneath. The 
electrons that are pushed out spiral along magnetic field 
lines and experience less scattering by heavy ions. Even- 
tually, they lose energy mostly by synchrotron radiation. 
As the electrons gradually gain speed also in the direc- 
tion of the magnetic field lines, the polarization of the 
emitted radiation also changes. 

Of course, only a complete MHD simulation can give 
precise results, showing the geometry of currents respon- 
sible for such a magnetic field, and could predict pre- 
cisely how much energy could be transformed into radi- 
ation, what would be the spectrum and how the polariza- 
tion would change with time. The above arguments may 
only give us some confidence that such processes actu- 
ally can occur at the Galactic centre. 

A SIMPLE SCENARIO 

Let us propose the following scenario to explain flares 
from the Galactic centre. Let us consider an object 
(LMS) on a highly eccentric orbit with periastron below 
~ 10 Tg, that got there by non-resonant relaxation. The 
tidal force melts the LMS and eventually brings it on a 
^ ~ 1 critical orbit. Just for the sake of simplicity, let us 
assume that the critical orbit is a parabolic one, with the 
critical radius of 4 rg. Let us assume moreover that the 
energy released during this process comes from the grav- 
itational potential energy of the object and is thus propor- 
tional to its mass. Since potential energy differences on 
an orbit at r = 4 rg are of the order of a few percent of 
mc^, it follows that objects producing the flares proba- 
bly have masses of the order of 10^" g. If the sources 
of flares were gaseous blobs of such a large mass, they 
would find themselves below the Roche radius far away 
from the black hole, and would therefore be completely 
disrupted before producing any modulation of the light 
curve. We conclude that the source of the flares is a small 
and initiaUy soUd object that orbits the black hole above 
the effective Roche radius, and is then brought close to 
the black hole, where it melts. 

As the tidal deformations on this orbit grow exponen- 
tially with time, the same happens with the magnetic field 
and the magnetic pressure, which resist the tidal squeez- 
ing and stretching - the magnetic field helps prevent- 
ing the object from falling apart. The accelerations of 
the electrons are the highest at the endpoints of such a 
thread, because the magnetic field is highly inhomoge- 



neous there. Consequently, most of the synchrotron radi- 
ation is emitted from that region. 

Since the tidal deformation tensor is traceless and the 
magnetic pressure tensor is not, the tidal squeezing of the 
object is stopped by magnetic pressure, while the tidal 
stretching along the orbit continues. Eventually, due to 
this imbalance, the object breaks into a large number 
of smaller pieces, just as in case of Plateau-Rayleigh 
instability. Every new piece emits most of the radiation 
at its endpoints. These smaller pieces move on the same 
orbit, because the magnetic field prevents them from 
completely spreading around the black hole. 

The light-curve of such an event can be constructed in 
the following way: the small, radiation emitting pieces 
are represented by A'^ small rigid objects, orbiting on a 
critical orbit with the critical radius of 4 r„. The orbit 
is chosen in such a way, to match the time-scales with 
the observations, i.e. 50 rg/c. Each of the A'^ objects has 
the same light-curve with a different phase shift. The ob- 
jects which are farther away from the black hole are less 
deformed, have lower magnetic field and thus lower lu- 
minosity with respect to the ones closer to the black hole. 
Due to exponentially increasing tidal deformation, the 
light-curves of all the objects are modified by an expo- 
nential function, to reflect this. Finally, all the modified 
light-curves are integrated into one. 

The appearance and the luminosity of this event is 
calculated with respect to the observers at different in- 
clinations and position angles with respect to the or- 
bit of the in-falling object using ray tracing techniques. 
Our simple scenario has five parameters: the mass of the 
black hole M^h, the tidal "squeezing" time scale Th, the 
length of the tidal tail Ztaib the inclination of the orbit ; 
and the longitude of the fine of nodes Q.. An example 
of such a simple model is a fit to an observed Galac- 
tic flare light-curve (flare of April 4,2007, [53]), shown 
in Fig. 7. This light-curve is consistent with the follow- 
ing parameters (assuming a Schwarzschild black hole at 
GC): nihh = 2.9 x IO^Mq, Th = 4.8min, /tail = 6rg, inch- 
nation near 90° and Q. is undetermined, since there are 
only two peaks in the light curve. The mass following 
from the fit does not quite agree with the measured one. 
The discrepancy can easily be explained by introducing a 
Kerr black hole with angular momentum parameter 0.4. 
However we prefer not to dwell on such details before a 
more detailed understanding is reached. 

CONCLUSIONS 

Given for granted that our Galactic centre harbours a 
super-massive black hole, one of its puzzling charac- 
teristics are the flares detected at various wavelengths. 
Their main characteristics are: a) short duration; b) quasi- 
periodic oscillations with the period of 17-22 minutes; c) 
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FIGURE 7. XMM-Newton/EPIC flare of April 4, 2007, fit- 
ted with the above simple model. The inset shows (roughly) the 
intensity at a fixed position of the image as a function of time; 

the rising beginning corresponds to the "squeezing" time-scale 
and the duration of the falling edge corresponds to the length 
of the tidal tail. 

total energy release of the order 10^'^ erg; d) strong in- 
dications of the presence of magnetic fields up to 100 
gauss. 

The exact cause of such flares is still unclear, and sev- 
eral models have been suggested to explain their ori- 
gin. We have proposed and explored the idea that the 
flares are produced by the final accretion of dense objects 
(LMS) with mass of the order of 10^ g onto the massive 
black hole in SgrA*. To this aim, we have investigated 
the effects of black hole tides on small solid objects in 
the vicinity of a massive black hole. Such objects are ex- 
pected to populate the Galactic centre as a result of being 
stripped off their parent stars. In [27], we have shown that 
solid objects, in the mass range ^ 10^^ — 10^^ g, melt in 
the vicinity of the black hole. Their subsequent orbital 
evolution naturally leads to capture orbits, most likely 
via unstable circular orbits. In an exhaustive numerical 
simulation we studied the evolution of shapes and light- 
curves as the objects progress toward the horizon of the 
black hole. We find that tides generally elongate objects 
into long thin threads, which extend exponentially on a 
time scale Tc w 11.3 rg/ c, if the orbital criticality param- 
eter ^ is close to 1. Moreover,if the object is bigger than 
> 0.3 Tg, it is completely disrupted before producing any 
significant modulation of the light curve. Therefore the 
object should be relatively small to produce the observed 
modulation. 



We studied numericaUy the tidal evolution of the 
shapes and the images of such objects in the gravita- 
tional field of a Schwarzschild black hole. Because of 
the complexity of the problem, we did not consider the 
more general Kerr space-time and limited the analysis to 
pressure-less gas. We find, quite generally, that such ob- 
jects are eventually elongated into thin threads, which on 
^ ~ 1 orbits extend exponentially. 

We point out that during such exponential elongation, 
conditions for magnetic flux conservation are met, so that 
relatively large magnetic fields can be generated. The ex- 
ponential growth of magnetic field density also allows 
the betatron mechanism which can generate highly rela- 
tivistic electrons. We propose this process as a possible 
source of galactic flares. We do not have the exact mecha- 
nism to extract radiation from gravitational potential en- 
ergy, yet we are able to model fight curves of galactic 
flares with simple assumptions. 

Finally, we would like to point out that our scenario 
differs from other models in the following respects: 1) 
Flares originate from marginally bound and not from 
marginally stable orbits. This brings them closer to the 
black hole and provides a mechanism for fast extraction 
of energy. It also agrees with recent measurements of the 
size of flaring region in the Galactic centre, showing that 
it is not larger than 6 rg, see e.g. [48], and relaxes con- 
straints on the angular momentum of the black hole. 2) 
Our model does not require a fairly rapidly rotating black 
hole as needed e.g. by [8, 49, 16]. 3) An accretion disk at 
the Galactic centre is not needed. In [50] it is highlighted 
that "there is no definitive evidence that Sgr A* has a 
disk". 4) Elongated structures invoked in other models 
to fit flare data, e.g. [51, 25, 52], develop naturally in 
the proposed model as a result of the tidal evolution of 
a melted body. 5) The evolving elongated structures also 
provide a natural mechanism for fast changes in mag- 
netic field density inferred from observations, increasing 
first as the body is squeezed by tides and finally decreas- 
ing as the body is crossing the horizon of the black hole. 
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